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SUN  workstation,  seismic  data  processing 


A  SUN-l/lOOU  graphic  workstation^i^pplied  by 'DARPA^ias  bee^  installed  in 
the  Selsmolc'gy  Computing  Lab^ratery>^at  Lamont-Doherty,  ana^has  become  the 
centerpiece  of  a  digital  waveform  processing  and  analysis  facility.  A  local 
urea  network  (LAN)  comprising  Ethernet /TCP-IP  protocols  connects  the  SUN-1 
with  the  Observatory's  pre-existing  VAX  11/780  and  PDP  11/70  Installations. 

The  operating  systems  on  the  VAX  (WS)  and  the  PDP  (Unix  bsd  2.9)  were  updated 
and  modified  to  provide  support  for  the  full  set  of  Ethernet  utilities.  The 
Ethernet  LAN  is  undergoing  start-up  testing  and  soon  will  be  fully  operational. 

While  the  SUN-1  offers  only  rudimentary  computational  po^er,  its  stare-of- 
the-art  graphic  capabilities  have  enabled  us  to  develop  a  ’'^ick-look*'algital 
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^data-handli ng  facility.-  Once  the  network  is  fully  operational,  the  resources 
of  the  Observatory's  other  computer  laboratories  will  become  available  to  the 
SUN-1.  The  VAX  will  be  used  as  the  network  flleserver,  peripheral  handler, 
and  computational  node,  while  the  SUN  will  be  used  for  the  graphic  analysis 
of  waveforms  and  digital  images.  Preliminary  experiments  with  the  network 
suggest  that  powerful  graphic  workstations  Integrated  vertically  with  file- 
servers,  computational  nodes,  and  supercomputer  communication  capability 
provide  a  flexible,  powerful,  and  productive  computational  environment  for  - 
research  in  seismology  and  tectonophysics.  IffuiDo S<-'/stn'ic.  c^ayo.., 
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IMTRODDCTION 

The  primary  purpose  of  this  contract  was  to  integrate  a  DARPA- 
supplied  Sun  workstation  with  the  existing  computer  facilities  avail¬ 
able  to  the  Seismology  group  at  Lamont-Doherty.  This  has  been  the 
first  step  in  greatly  expanding  Lamont-Doherty ' s  capabilities  for  the 
collection,  storage,  retrieval,  display  and  analysis  of  digital 
seismic  waveform  data.  The  work  covered  in  this  contract  primarily 
involved  hardware  installation  and  modifications /add it  ions  to  oper¬ 
ating  systems  to  allow  direct  communication  between  computers. 

Some  software  development  has  been  done,  but  this  activity  is 
primarily  taking  place  under  a  new  DARPA  contract.  The  main  software 
oriented  goals  are  the  development  of  interactive  graphic  software  for 
the  display  and  preliminary  analysis  of  digital  seismic  waveform  data 
and  the  development  of  interactive  modeling  capabilities  for  more 

sophisticated  waveform  analysis. 

The  funds  under  the  contract  provided  for  the  personnel  time  for 

the  installation  of  the  Sun  system  and  purchase  of  limited  hardware 

for  an  Ethernet  interface  between  the  Sun  and  existing  computer 

systems . 


NBTHORKIIIG 

A  block  diagram  showing  the  existing  and  proposed  developments  in 
the  computing  network  at  Lamont-Doherty  is  shown  in  Figure  1.  Indi¬ 
cated  on  this  figure  are  those  components  that  have  resulted  directly 
from  this  contract  and  those  that  have  been  provided  by  the  Obser¬ 
vatory,  primarily  as  a  result  of  stimulation  from  this  work. 

The  workstation  as  provided  by  DARPA  consisted  of  a  SUN-l/lOOU 
with  an  80  Mbyte  disk  and  a  1600  bpi  tape  drive.  The  system  had  a 
number  of  limitations  which  initially  made  it  difficult  to  realize  its 
full  potential  as  a  seismic  analysis  workstation.  The  most  critical 
limitations  were  the  lack  of  a  local  printer/plotter,  limited  disk 
storage,  relatively  slow  CPU  (no  arithmetic  processor)  and  differences 
in  graphics  primitives  and  screen  size  from  current  Sun  models. 

The  initial  step  to  overcome  these  limitations  was  to  link  the 
Sun  to  the  PDP  11/70  via  uucp  at  9600  baud.  This  provided  the  Sun 
with  access  to  the  PDP's  two  Versatec  printer/plotters,  its  disk  stor¬ 
age,  and  its  connections  with  the  outside  world.  The  uucp  link  has 
been  most  used  for  printing  and  plotting  as  there  is  no  local  printer 
on  the  Sun.  Uucp's  remote  execution  capabilities  have  allowed  the  use 
of  a  simple  command  on  the  Sun  (vpr  filename)  to  combine  the  file 
transfer  and  printing  operations.  The  same  is  true  of  Sun  graphic 
screen  dumps.  The  raster  files  on  the  Sun  are  transferred  to  the  PDP 
for  plotting  on  the  Versatecs.  With  uucp  providing  the  transfer  and 
remote  execution,  a  program  to  actually  write  the  raster  file  to  the 
plotter  was  written  locally. 

While  our  original  intention  was  to  establish  an  Ethernet  link 
between  the  Sun  and  our  Seismology  PDP  11/70,  we  have  been  able  to  use 
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this  development  as  leverage  to  convince  the  Observatory  to  provide 
the  funds  to  begin  the  installation  of  a  campus-wide  Ethernet  link, 
starting  with  a  link  bettieen  the  Observatory's  main  VAX  11/780  and  the 
Seisaiology  computer  facility.  The  Observatory's  VAX  operates  under 
the  DEC  WS  operating  system.  To  facilitate  communication  with  the 
UNIX-based  11/70  and  Sun  computers  in  Seismology,  a  UNIX  emulator, 
EUNICE,  has  been  installed  (with  Observatory  funds)  on  the  VAX,  pro¬ 
viding  full  UUCP  and  Ethernet  protocols.  A  local  area  Ethernet,  a 
major  segment  of  it  using  fiber  optic  cable  and  again  funded  primarily 
with  Observatory  funds,  will  connect  the  Sun,  the  PDP  11/70  and  the 
VAX  11/780.  The  VAX  archives  Lamont-Doherty  Geological  Observatory's 
(logo's)  geophysical  data  bases  and  supports  extensive  peripherals 
such  as  image  processing  hardware,  laser  printing,  and  large  pen 
plotters.  All  these  resources  will  become  available  to  the  Sun  over 
the  nettrork.  The  wiring  for  the  connection  between  the  Sun  and  the 
PDP  has  been  installed  and  the  software  on  the  PDP  is  now  in  the 
process  of  being  upgraded  to  allow  networking.  (We  acknowledge  sub¬ 
stantial  assistance  from  Keith  Bostic  at  CSS  and  Scott  Bradner  at 
Harvard.)  The  VAX-1 1/70  link  is  now  operating  through  uucp  at  9600 
baud.  The  VAX  will  join  the  full  Ethernet  network  as  soon  as  fiber 
optic  repeaters  arrive. 

The  network  will  make  the  CPU  power  of  the  VAX  easily  available 
to  the  Seismology  computing  facility,  alleviating  the  CPU  limitations 
of  the  Sun-1.  The  disk  storage  available  to  the  Sun  is  also  being 
improved  by  the  purchase,  again  with  Observatory  funds,  of  one  400 
Mbyte  disk  drive.  The  next  step  in  expanding  the  network  is  an  addi¬ 
tional  Sun  color  workstation  and  a  proposal  has  been  submitted  to  NSF 
for  this  purpose. 

SOFTWARE  DKVELOFHEVr 

We  are  developing  an  interactive  seismic  analysis  system  which  can 
be  implemented  on  a  local  area  network  of  high-powered  graphic  work 
stations.  The  SUN-1 /lOOU  on  loan  from  DARPA  has  been  the  cornerstone 
of  this  emerging  system,  as  well  as  a  test-bed  for  software  develop¬ 
ment.  In  the  course  of  this  effort,  we  have  identified  four  cate¬ 
gories  of  software  development  tasks: 

1.  Interactive  catalog  searches  and  retrieval  of  GDSN  waveform 
data:  As  a  designated  regional  data  center,  LDGO  is  a  repository  for 
ISC  and  PDE  catalog  data,  ISC  parameter  data,  and  GDSN  network-event 
data.  We  have  developed  search  algorithms  for  these  data  bases,  but 
have  been  unable  to  keep  most  of  the  catalog  information  on-line 
because  of  disk  space  limitations.  In  addition,  we  have  developed 
interactive  routines  for  the  playback  of  GDSN  network-event  tapes 
(Figure  2).  We  are  developing  management  software  to  integrate  the 
ISC  and  PDE  catalog  data  with  a  catalog  of  the  events  contained  in  the 
GDSN  data  base,  so  that  search  procedures  can  be  integrated  with  the 
waveform  data  playback. 

2.  Interactive  waveform  parameter  extraction,  and  interactive 
parameter  and  waveform  modelling  for  studies  of  source  and  structure: 
We  have  installed  interactive  parameter  extraction  software  (spex). 
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provided  by  Science  Horizons  (Wang  et  al . ,  1983).  In  addition,  we 
have  rewritten  an  interactive  velocity  modelling  code  (originally 
written  by  P.  Shearer  of  Scripps  Institution  of  Oceanography)  to  con¬ 
form  to  the  SIGGRAPH  Core  standard,  and  have  installed  it  on  the  SUN 
(Figure  3).  This  combination  of  codes  gives  us  the  capability,  for 
example,  to  forward-model  travel-time  parameters  extracted  from  wave¬ 
form  data  bases  quickly  and  nearly  automatically.  We  will  expand  the 
capabilities  of  the  velocity-modelling  code  to  include  preliminary 
travel-time  inversion  and  waveform  modelling  capabilities.  The 
graphics  capabilities  of  the  SUN  are  instrumental  to  the  development 
of  an  efficient  user  interface  to  these  codes. 

3.  Synthetic  seismogram  codes  and  software  development  systems: 
While  the  SUN-1  does  not  possess  a  processor  with  sufficient  power  to 
calculate  synthetic  seismograms  for  realistic  earth  models,  its 
multi-window  environment  and  UNIX  operating  system  provide  for  excep¬ 
tionally  rapid  development  and  debugging.  Once  debugged,  complex 
codes  are  uploaded  to  a  faster  processor  for  computation.  Results  are 
downloaded  to  the  SUN  for  graphic  analysis  (Figure  4).  We  are  imple¬ 
menting  an  organizational  framework  for  software  development  that  will 
permit  us  to  keep  track  of  the  multitude  of  codes  and  revisions  pro¬ 
duced  in  a  research  environment. 

4.  Interactive  graphic  routines  for  plotting  2-D  and  3-D  geo¬ 
physical  data  sets:  We  have  written  and  installed  on  the  SUN  a  set  of 
command- format  plotting  codes  for  2-D  and  3-D  gridded  geophysical 
data.  These  currently  are  being  debugged.  We  have  also  developed  the 
capability,  using  SUN  graphics  primitives,  to  produce  publication- 
quality  line  and  gray-shade  plots  of  3-£)  geophysical  data  (Figure  4). 
We  are  evaluating  the  possibilities  of  translating  these  codes  into 
the  CORE  graphics  standard. 

Some  of  these  development  activities  might  seem  to  reproduce  the 
functions  of  CSS,  but  we  are  confident  that  the  software  developed  at 
L-DGO  will  be  of  interest  to  CSS.  While  we  recognize  the  importance 
of  CSS  as  a  national  data  center,  we  believe  that  the  development  of 
^  local  data  base  management  capabilities  and  analytical  software  will 

I  complement,  rather  than  compete  with,  facilities  offered  by  CSS.  Not 

only  will  our  local  capabilities  enhance  the  productivity  of  our 
investigators  doing  research  of  interest  to  the  DARPA  program  for 
treaty  verification,  but  we  envisage  that  the  development  of  local 
capability  integrated  vertically  with  CSS  communication  capability 
will  expand  the  usage  of  CSS  facilities  by  our  scientists. 

ADDITIOHAL  ACTIVITIES  BSUITBD  TO  THIS  COHTRACT 

Digital  topographic  data  for  Soviet  Central  Asia  were  requested 
under  this  contract  for  use  in  comparing  topography  and  seismicity. 
These  data  were  supplied  during  May-July  1985.  Two  hundred  cells  of 
l*xl*  each  with  1200  points/degree  (appox.  100  m  spacing;  total  of 
3x10®  points)  are  included  in  the  10*x20*  area  covering  the  Tien  Shan, 
Pamir,  Hindu  Kush  and  Tadjik  Depression.  These  data  are  now  being 
entered  into  a  format  suitable  for  comparison  with  the  available 
seismicity  data  and  for  analysis  with  our  image  processing  facility. 


-  u  *  «  ^ 


The  Palisades  seismograph  station  (PAL)  at  Lamont -Doherty  has 
recently  been  upgraded  to  digital  recording.  The  data  are  collected 
on  a  dedicated  multibus  microprocessor  and  written  continuously  to  9 
track  tape.  The  tape  is  first  processed  on  the  Sun  and  then,  as 
necessary,  events  are  transferred  to  the  PDP  or  VAX  for  further  pro¬ 
cessing.  A  program  has  been  written  to  decode  the  9-track  tape,  which 
is  in  logo's  own  format,  into  CSS  'wfdisk'  format  and  selectively 
retrieve  the  events.  A  second  program  selectively  extracts  events  to 
be  saved  on  an  archive  tape. 

RBFKRKIKES 

Wang,  J.,  R.  M.  Rother,  V.  L.  Hutchison,  J.  H.  Alexander,  and  J.  B. 
Minster,  Seismic  parameter  extraction  package,  S-Cubed  Report 
#SSS-R-84-6423,  1983. 


CHARGBS  TO  5-22107  (NETWORK  RELATED) 


1  Fiber  optic  cable  for  Ethernet  repeater  $1750.00 

2  Interlan  Ethernet  transceivers  580.00 

2  DEC  Ethernet  transceivers  600.00 

4  Transceiver  cables  304.00 

4  Ethernet  terminators  75.00 

Tools  253.58 


TOTAL 


$3562.58 
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